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Abstract. 

As part of the HSF/SFIS GRB host survey program we present the detection of three faint gamma-ray burst 
(GRB) host galaxies based on an accurate localisation using ground-based data of the optical afterglows (OAs). A 
common property of these three hosts is their extreme faintness. Fhe location at which GRBs occur with respect 
to their host galaxies and surrounding environments are robust indicators of the nature of GRB progenitors. 
Fhe bursts studied here are among the fou r most extreme ou tliers, in terms of relative distance from the host 
center, in the recent comprehensive study of Bloom et al. ( :2002| ). We obtain a revised and much high er probability 



that t he galaxies identified as hosts indeed are related to the GRBs (P(n c hanco)= 0.78, following Bloom et al 



2002)), thereby strengthening the conclusion that GRBs are preferentially located in star- forming regions in 
their hosts. Apart from being faint, the three hosts consist of multiple structures, indicative of merging and active 
star-formation. For one of the hosts, GRB 980329, we estimate a photometric redshift of z ~ 3.5. 



Key words. Gamma rays: bursts - Cosmology: observa- 
tions - Galaxies: starburst 



* Based on observations made with the NASA/ESA Hubble 
Space Felescope, obtained from the data archive at the Space 
Felescope Institute. SFScI is operated by the association of 



1. Introduction 

There are two, possibly fundamentally different, classes 
of gamma-ray bursts (GRBs) - short and long duration 
bursts (being shorter and longer than ~ 1 second respec- 
tively). So far every GRB with an identified optical af- 
terglow (OA) belongs to the class of long duration bursts, 



Univer: 

contrac t NAS 5 26555. Baaed on obacr - vationa made with ESQ 



with the possible exception of GRB 000301C (Jensen et al 



ities fui Research in AsQuuuiuy, Inc. uudei the NA3A — onnl K mi . i x r j i c j — 

,r . ~ ^ . , ™^ 2001). 1 here are two general types of models for produc 



Telescopes at the La Silla or Paranal Observatories under pro- 
gramme ID 60.A-0552. Based on observations made with the 
Nordic Optical Telescope, operated on the island of La Palma 



ing GRBs, the first involves the merging of binary com- 



pact stars (Paczyhski 1986), the second is related to the 
death of very massive stars (Woosley 1993; MacFadyen & 



jointly bv Denmark, Finland, Iceland, Norway and Sweden, in 
the Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astrohscia de Canarias 



Woosley 1999; Vietri fc Stella 1998). The observations of 



SN bumps, Fe K-line and OA localizations close to star- 
forming regions favour the latter type of models. Due to 
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the relatively short lifetime of such massive stars, one 
expects them to exist in or very close to star-forming 
regions where they are born. Thanks to the extensive 
observational efforts in pursuing GRB events at various 
wavelengths there are now several examples where such 
positional correlations between the GRB OA and star 



forming regions in the host are seen, for instance Fynbo 
et al. pOOOj ) ( GRB 980425), |Holland fc Hjorth (1999^ 



(GRB 99012 3 ), [Hjorth et al. (2002| ) (G RB 980613), [Bloom 



et al. ( |200lD , |Chary et al. (2002[ ) and [Frail et al. (2002| ) 



(GRB 010222). There are, however, a few cases in which 
the GRB does not seem to originate from intense star- 
formin g reg ions (SFRs), eg. and GRB 990705 (Andersen 
et al. 2ft02[ ). 

In this paper we present localisations and host 
candidates of three OAs based on ground data and 
Hubble Space Telescope (HST) STIS imaging data 
from the Cycle 9 program GO-8640 "A Survey of 
the Host Galaxies of Gamma-Ray Bursts" 



( Holland 

et al. 3000| ) (data and further information available at 



http : //www . if a. au. dk/~h.st/grb_hosts/index .html). 

In Section ^ we describe the image processing applied to 
the data, in Section || the astrometry resulting in the OA 
STIS-image localisations is described. In Sections [| - || 
we indentify the hosts, present photometry and discuss 
the host environment. Finally, in Section |?] we discuss the 
implications of our localisations and host identifications. 
Specifically, we re-compute using our host identifications, 
the probability computed in Bloom et al. (2002) that 



none of the host identifications in that sample are random 
galaxies. Finally, we summarize our own results. 



2. Data 

Using the STIS-instrument onboard the HST the GRB- 
systems were observed using a four-point dithering pat- 
tern with shifts of 2.5 pixels (~ 0'/127) between exposures. 
The data was pre-processed using the standard STIS 
pipeline and combined using the DITHER (v2.0) soft- 
ware (Fruchter & Hook 2002) as implemented in IRAFQ 
(v2.11.3) and STSDAS (v2.3). The STIS images were driz- 
zled using 'pixfrac=0.6' and 'scale=0.5' (giving a pixel size 
of 0'.'0254). Note that drizzling introduces correlated noise 
between neighbouring pixels. All GRB-systems were ob- 
served using the STIS 50CCD (hereafter CL) passband 
with pivotal wavelength PivW= 5851.5A. GRB 980329 
was also observed using the STIS F28X50LP (hereafter 
LP) passband, PivW= 7228.5 A. 

The photometry was performed in circular or elliptical 
apertures, as appropriate, according to the morphology of 
the host. The size of the apertures was selected so as to 
measure the total flux, by first choosing a plausible shape 
and then increasing the size until no gain in flux could 
be achieved. The sky was measured in an annulus with 



1 Image Reduction and Analysis Facility (IRAF), a soft- 
ware system distributed by the National Optical Astronomy 
Observatories (NOAO). 



corresponding shape to the object aperture and with inner 
and outer annulus 1.5 and 4 times the object aperture. 
For the STIS zero-points, we adopted the values found by 
pardner et al. (2000j ) for the HDF-south. The zero-points 
used were therefore ZP CL = 26.387 and ZP LP = 25.291. 
Foreground (Galaxy) extinction estimates were computed 
using the on-line NED extinction calculator^ based on the 
dust maps provided by [Bchlcgcl et al. (199^ ) . 

Signal-to- noise (S/N) estimates of the host detections 
were computed as the ratio between the measured counts 
in a circular aperture centered on the object and the sky 
variance as measured from 10 similar apertures at random 
positions (on the sky). The (circular) aperture diameter 
used was 19, 19, and 9 for GRB 980329, GRB 980519 and 
GRB 990308, respectively. Note that the errors on the pho- 
tometry is computed as in IRAF DAOPHOT/APPHOT 
and does not necessarily correspond to the S/N estimates. 

The host pixel positions were determined by using the 
IRAF/APPHOT CENTER task and the 'centroid' algo- 
rithm therein. The 'centroid' algorithm computes the in- 
tensity weighted means of the marginal profiles in x and 
y. The results are given in Table [l| 

Throughout this paper we use the following cosmolog- 
ical parameters; = 0.3, Qa = 0.7, H = 70 km s _1 
Mpc" 1 . 



3. Astrometry 

The procedure used for relative astrometry from ground- 
based images of the afterglow to the STIS data depends 
on the available ground-based data, but for all three ob- 
jects least squares affme transformations were used. An 
affine transformation is the simplest possible transforma- 
tion which allows for deviations from square pixels in ref- 
erence and target images. It is known that some of our 
ground-based reference data were obtained with CCDs 
which have slightly non-square pixels. As we were not able 
to establish any clear correlation between position and 
astrometric residuals, it is justified to keep all transfor- 
mations strictly linear, as is also preferred with a limited 
number of tie objects. The accuracy by which the position 
of a (point) source in a well sampled image can be deter- 
mined is dependent on the signal-to- noise ratio (S/N) with 
which the source is detected, and the full width at half 
maximum (FWHM) of the point spread function (PSF). 
For low S/N the astrometric standard error per axis can 
be approximated with cr pos = <jpsf/(S/N), where ctpsf is 
the standard deviation of the Gaussian approximating the 
PSF, which formally equals FWHM/2.35. For high S/N, 
accuracy is limited by errors in the detector pixel geom- 
etry, usually at the level of 1/20 of a pixel or less. If the 
S/N is significantly above 20 or the FWHM is sampled 
by less than about 3 pixels, this approximation of the as- 
trometric error is not valid. Whenever several individual 
images of the afterglow is available, a transformation to 
the STIS reference frame is established for each image. 



http:/ /nedwww. ipac.caltech.edu/ 
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Table 1. GRB R.A. and Dec. and the best fit pixel localisation (incl. 1-a errors) of the OA in the CL images. For 
each OA, three reference objects are given including their pixel positions 



Target 


R.A. Dec. 


OA 

X Y 


Host gala 
R.A. Dec. 


X Y 


980329 
ref. f 
ref. 2 
ref. 3 


07:02:38.02 +38:50:44.3 
07:02:38.f5 +38:51:03.7 
07:02:37.43 +38:50:34.2 
07:02:38.97 +38:50:33.2 


999.36 ±1.34 1032.33 ±1.22 
941.32 1800.06 
1274.41 635.89 
564.13 598.21 


07:02:38.07 ±38:50:44.3 


995.55 1027.35 


9805f9 
ref. 1 
ref. 2 
ref. 3 


23:22:21.54 +77:15:43.2 
23:22:17.78 +77:15:51.3 
23:22:20.25 +77:15:23.8 
23:22:27.36 +77:15:43.7 


1039.19 ±0.68 990.77 ±0.71 
416.04 1128.28 
1410.45 1673.18 
1514.40 462.34 


23:22:21.34 ±77:15:43.7 


1045.27 975.98 


990308 
ref. 1 
ref. 2 
ref. 3 


12:23:11.49 +06:44:04.7 
12:23:10.06 +06:43:50.0 
12:23:12.30 +06:44:24.3 
12:23:10.95 +06:44:20.8 


813.46 ±2.79 922.50 ±2.03 
1651.271 341.10 
340.29 1698.05 
1131.42 1558.46 


12:23:11.49 ±06:44:04.7 


811.70 921.80 



As the error in the transformed afterglow position in the 
STIS image is always completely dominated by the as- 
trometric error in the ground-based image, the errors of 
individual transformed afterglow positions are in practice 
independent. The standard deviation of the transformed 
positions is therefore an estimate of the actual astrometric 
error. 



4. GRB 980329 

GRB 980329 was detected by the BeppoSAX satellite on 
1998 March 29.16 UT. The radio and optical counterparts 



were discovered by Taylor et al. (1998) and Djorgovski 



ct al. (1998), respectively. The latter claimed the detection 
of the host galaxy at an apparent magnitude of R ~ 25.7 



Palazzi et al. (1998[) jGorosabel et al. (1999] ) , |Reichart et al 



(1999| ) presented optical and near infrared detections of 



the OA and found it to decay at a rate typical for most de- 
tected OAs (ao ~ 1-2). This decay slope is in good agree- 
ment with that found in X-rays, where ax = 1-35 ± 0.03 
( in 't Zand et al. 1998 ). The maximum measured (early) 
brightness of the OA was I = 20.8 and R= 23.6 , leading 
to an extre mely red colour R— I = 2.8 ± 0.4 ( Reichart 

(NIR) observations 



et al. 



199E ). Near- infrared (NIR) observations, on the 
other hand, showed that the NIR c olours are appr oxi- 
matcly flat. These measurements led Fruchter (1999 ) to 



argue that the red colour could be caused by the Ly-a for- 
est if the redshift was z ~ 5. This claim was subsequently 
challenged by the apparent non-detection of a Ly-a forest 



in Keck II spectra (Djorgovski et al. 2001). Recently, Yost 
et al. (2002) presented supplementary multi- wavelength 



broad-band photometry of this burst and claim to rule out 
z > 5 based on an afterglow model fit to the data. The 
NIR photometry shows that the HST/NICMOS October 
1998 data (GO-7863, PI: A. Fruchter) contains a signif- 



We retrieved NTT/EMMI R-band images of 
GRB 980329 from the ESQ archive, obt ained on 
March 29.99 and 30.99 QPalazzi et al. 1998j ). As the 
afterglow was detected at low signal-to-noise in the 
late-time images, the astrometry was in this case derived 
from the combined image. Seven tie objects were used 
for the astrometric solution. The astrometric error, 
as estimated from the residual of the tie object fit, is 
estimated to be about 1.25 drizzled STIS pixels, or 0'.'03, 
which should be compared with an expected error of 
0'.'02, as estimated from the S/N of the OA image. The 
error in the transformation from STIS CL to the STIS 
LP image is a small fraction of a drizzled STIS pixel and 
can therefore be ignored. The best fit localisation in the 
CL image is given in Table 

An excerpt of the STIS CL and LP images are shown 
in Fig|l] centered on the host and indicating the OA posi- 
tion. In the CL image several unresolved knots are seen on 
top of a low surface brightness area within an aperture of 
0'.'5. The measured ABMAG in the CL-band within this 
aperture is 27.5 ± 0.2. In the LP image an extended ob- 
ject is seen, but the knots seen in the CL-band are not 
detected. We find 26.6 ± 0.2 within the same aperture in 
the LP-band. Foreground extinction corrected photometry 
and detection signficance estimates are given in Table f§ 
Photometry of the three brightest knots yields a total 
magnitude of 28.1 ± 0.1 in the CL-band and 27.7 ± 0.2 
in LP. The flux in the LP measurement stems primarily 
from the underlying galaxy complex. Within an arc-second 
of the OA position at least two fainter extended objects 
or structures are seen to the North and North-East in the 
CL-band. Their distance relative to the OA position is 
approximately 1 and 0.6 arc-second and photometry mea- 
surements yield 28.5 ± 0.2 and 29.5 ± 0.5, respectively. 

The STIS CL and LP measurements in addition to 
the Keck/ESI R = 26.53 ± 0.22, I = 26.28 ± 0.27 ([Bloom 



icant c ontribution from the OA and therefore it is not et al. 2002[ [Yost ct al. 2002|) and NIRC K = 23.04 ± 0.42 



suited for host photometry. The GRB has also been found 



to be heavily extincted by dust (see Lamb et al. 1999; 



Bloom et al. 2002; Yost et al. 2002). 



measurements (Yost et al. 2002) provide an excellent op- 
portunity to estimate the photometric redshift of the host. 
Using the Bayesian photometric redshift (BPZ) estimation 
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Fig. 1. Sub-section of the STIS CL and LP images centered on GRB 980329. A dashed ellipse indicates the 3-u OA 
localisation error while the solid circle gives the aperture used for photometry and a cross in the CL image marks the 
host center. The dash-dotted ellipse marks the 3-cr error circle of Bloom et al. All images were smoothed by a 1.1 pixel 
Gaussian using IRAF/GAUSS. Axis units are in arc-seconds, North is up and East is to the left. 



software of Bcnitcz (2000| ) and restricting z > 1 (due to ab- 
sence of expected emission lines in spectra of the host, see 
Yost et al. 2002] ) we find z ~ 3.6 with the best fitting SED 
corresponding to an Im galaxy type. Redshifts of z < 1.2 
and z > 4.2 are excluded at the 95% confidence level (z 
> 5 excluded at 99.99% level). These redshift estimates are 
consistent with the constraint z < 3.9 based on the non- 
detection of the Ly-a forest in a Keck II spectrum of the 



host galaxy (Djorgovski et al. 2001) and a far- ultraviolet 
extinction curve constraint giving 3 < z < 5 (Lamb et al 

It is puzzling that the CL-band shows a clear multi- 
component nature, with at least three unresolved knots 

within 0.5", whereas the LP-band docs not. We find an We used the NOT observations of the OA fljaunscn 



( Jaunsen et al. 19"98| ) . Jaunsen et al. (2001 ) presented a 
homogeneous optical data sample obtained at the Nordic 
Optical Telescope (NOT) and found that the light-curve 
breaks at around to + 0.5 days, with a pre- and post-break 
power law of aoi = 1-73 and ao2 = 2.22, respectively. The 
X-ray (pre-break) power law was approximately, ax = 1.6 
though this estimate is somewhat uncertain as it is based 
on combining the BeppoSAX WFC and NFI measure- 
ments. The maximum measured flux of the OA was 1=18.4 
and R=19.8 at to + 0.346 and to + 0.536 days, respectively. 
Correcting for the decline, the approximate R— I colour 
was ~ 0.6. 



upper |imit colour for the knots of CL-LP < 0.25, while 
for the integrated colour of the host complex, we measure 
CL-LP w 0.8. The compact knots does therefore appear 
very blue, implying that Ly-a may be present in the CL- 
band. Our best estimate of the redshift is therefore zw 
3.5. T aking the m easured E 7 = 5.5 x 10~ 5 erg cm 2 flin 't| 
we find an isotropic gamma-ray energy 
10 54 erg. A ssuming a total av erage energy 
erg (Frail et al. 2001) we estimate 



Zand e t al. 1 
of E iso = 1.4 x 
for GRBs of 5 x 10 50 
a jet opening angle of 9q = 1.5 deg, indicating a highly 
collimated beam. 



5. GRB 980519 



GRB 980519 was detected by GRO/BATSE ( |Mullcr 1998Q 
on 1998 Ma y 19.51 UT and subsequent BeppoSAX WFI 
localisation ( Piro 199~8| ) enabled the detection of the OA 



et al. 2001 ) to derive an accurate position of the burst 
in the STIS image. As the GRB 980519 afterglow was 
discovered and followed up at very high airmass, differen- 
tial color refracti on (DCR) affects th e astrometry signifi- 
cantly. Following Monet et al. (1992 ) we correct for DCR 
by a term which is proportional to sec(Zdist) and depends 
on the color of each individual object. Since the colors of 
the tie objects and the afterglow cover a narrow range, the 
color dependence was approximated with a linear relation. 
The afterglow was observed in the R and the I bands, so 
two different DCR corrections must be applied. To mini- 
mize the number of free parameters, we choose to derive a 
theoretical estimate of the relative amplitude of the DCR 
correction in the R and I bands. From the different refrac- 
tion across the filter bandwidth it is found that DCR in 
R is about 2.65 times DCR in I. By using this factor, only 
the amplitude of DCR as function of color for the com- 
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Fig. 2. Sub-section of the STIS CL-band image centered on the OA for GRB 980519 and GRB 990308. Axis units 
are in arc- seconds, North is u p and East is to the left. Annotations as in Fig. |l| in addition to a dotted ellipse which 
marks the Bloom et al. C2002|) host identifications. 



plete data set of six images has to be determined. This 
is done by minimizing the standard deviation of the six 
independent localisations. For each localisation between 7 
and 10 tie objects were used. The difference between the 
final averaged positions with and without DCR correction 
is of the order of 0.7 drizzled STIS pixels, comparable to 
the l-cr astrometric error. In effect, by applying the DCR 
correction, the average of the three R-band positions co- 
incide with the average of the three I-band positions. The 
OA localisation result is given in Table [p. 

At the OA position we detect two extended ob- 
jects, clearly visible in Fig. |[ The photometry yields CL 
ABMAGs of 27.1 ± 0.1 and 27.9 ± 0.1 for the southern 
and northern component, respectively. The OA is located 
in the very outskirts of the northern component (S/N~9), 
where a faint blob coinciding with the position of the OA 
is detecte d (see also Table p|) . A host detection was re- 
ported by [Bokolov et al. (199q ) and Bloom et al. (199S ) at 



an estimated Cousins R VEGAMAG of 26.1 ± 0.3. This 
detection, however, consisted of the smeared sum of both 
objects. By using a larger aperture of 1'.'88 (enclosing both 
objects) we find ~ 26.45 ±0.10. Assuming a flat spectrum 
and using the STSDAS SYNPHOT/CALCPHOT we con- 
verted the CL ABMAGs to Cousins R VEGAMAG, giving 
~ 26.0, in agreement with the earlier combined detections. 
It is also worth noting that the combined flux of the two 
objects amounts to 26.7 (ABMAG) as compared to 26.45 
for the large aperture. The difference in flux can be at- 
tributed to the very low surface brightness (> 29) regions 
in the vicinity of the two major components. Assuming a 
rcdshift larger than 0.5, which is reasonable given the red- 
shift distribution of other GRB hosts and the faintncss of 



the host, we note that the average angular scale of 1 arc- 
second is ~ 7 ± 2 kpc for our assumed cosmology. Given 
this angular scale and the disk-like morphology it is most 
likely that the two detected objects are galaxies which are 
in the process of merging. This is supported by the low 
(< 0.003) integrated probability ( |Gardner et al. 2000] ) of 
having two objects of this brightness within ~ 2 sq. arc- 
seconds. The faint neighboring patches are therefore likely 
to be smaller galaxy fragments belonging to the merging 
system. 

6. GRB 990308 

GRB 990308 was detected by GRO/BATSE on 1999 
March 8.21 UT, on the RXTE All-Sky monitor and also 
weakly by the Ulysses GRB detector. An OA was detected 
using the QU EST camera on a 1.0 m Schmidt telescope 
in Venezuela (Schaefer et al. 1999). These optical mea- 
surements give a power law of a = 2 ± 2 (Schaefer et al 



1999), and are the only data in which the OA is detected. 



Early non-detections by LOTIS and Super-LOTIS sug- 
gest a < 1.3, while later non-detections by the WIYN 
and Keck telescopes set the constraint a > 1.1. Taking all 



constraints into consideration, Schaefer et al. (1999) found 
a best fitting constant power-law, a = 1.2 ± 0.1. 



We used the original QUEST data (Schaefer et al 



1999| ) to transform the OT position to the STIS clear 
image coordinate system via an intermediate transforma- 
tion to a combined NOT image due to the lack of com- 
mon tie objects between the QUEST and STIS images. 
The NOT/ALFOSC image was combined using standard 
tools based on 7 R-band images obtained specifically for 
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Table 2. Aperture photometry of GRB hosts in AB magnitudes and corrected for foreground extinction according to 
Schlcgel et al. (199^ ) using the Johnson V-band estimate for CL and Cousins I-band estimate for LP. 



GRB 


Obs. Date 


E(B-V) 


CL 


LP 




absolute 


days 




int. 


mag 


s/n 


int. mag s/n 


980329 


24/26 Aug 2000 


~ 880 


0.073 


8072 


27.2 ±0.1 


10 


8156 26.2 ±0.1 9 


980519 


7 Jun 2000 


~ 750 


0.240 


8983 


27.0 ±0.2 


9 




990308 


19 Jun 2000 


~468 


0.023 


7842 


29.7 ±0.4 


5 





this purpose on 29-30 March 2001 with a total integra- 
tion time of 6000s. The final STIS position is the average 
of the three individually transformed positions, based on 
the three useful QUEST images. The error is estimated 
from the scatter of the individual positions and includ- 
ing a contribution from the NOT to STIS transformation, 
with the latter being small. In total 11 tie objects were 
used in the transformation to the STIS image. The final 
error is consistent with the signal-to-noise of the OT detec- 
tions in the QUEST images and the equation introduced 
in Sec. ^. The final STIS position and error estimate is 
given in Table [lj At the locus of the OA we marginally 
detect a very faint point-like object (see Fig. ||) which we 
measure to have a STIS CL magnitude of 30.1 ± 0.4 with 
a detection significance of S/N ~ 5 (foreground extinction 
corrected photometry given in Table ||) . Including the faint 
extended emission north of this object gives 29.9±0.4. We 
also detect an extended object 0'.'3 to the south with an 
estimated magnitude of 29.8 ± 0.4 (S/N ~ 3) and 0'.'8 to 
the North a much larger disk-like object with a magnitude 
of 27.7 ± 0.1. Using the galaxy counts of Gardner et al 



(2000[j]we find a relatively low probability (~ 0.02 within 
a radius of l"of the OA position) that the three objects 
are projected neighbors. 

Could the point-like component coincident with the 
OA location in fact be the OA itself? Assuming a constant 
power-law one can deduce the decay slope by interpolat- 
ing the brightness at the time of the first V-band QUEST 
observation and the measured brightness in the STIS ob- 
servations 468 days after the burst. This gives a power-law 
exponent of a ~ 1.35, consistent with the best estimate of 
the power-law slope, a ~ 1.2 (based on all available data). 
If correct, this would be the latest trace to date of an OA, 
at to + 468 days. Another possibility is that the point-like 
component is caused by some re-brightening mechanism, 
such as eg. dust echoing. 

In summary, wc identify the object coincident with the 
OA localisation as the possible remnant OA (point- like) or 
the host (extended). If the point-like object turns out to be 
non-variable and therefore not the OA, then it must have 
been fainter than ~ 30 mag. This implies that the late 
time decay slope must have been larger than 1.35. This 
scenario and the constraints from the early data could 
be explained by introducing a break in the light-curve. 
Specifically, an early a ~ 1.3 slope (as supported by the 
early LOTIS data) followed by a steeper slope fits this 
scenario well. A revisit of this field with HST±ACS is 
required to disentangle these ambiguities. 



7. Summary 

We have localised the three OAs to high precision in the 
STIS images and identify the host as the nearest detected 
object of the OA position. The GRB 980329 host galaxy 
redshift is estimated to be z ss 3.5. For all three candidate 
hosts we detect faint extended structures within a radius 
of ~ 1" (~ 7 proper kpc for z > 0.5). This scale is sim- 
ilar to that seen between tidally interacting and merging 
galaxies (eg. Borne et al. 2000 ). The hosts show signs of 
sub-structure (possibly star-forming and/or merging ele- 
ments); 3-4 blue knots in the GRB 980329 host, a knot in 
the northern edge of the GRB 980519 host and a point- 
source within the 3-cr localisation error of GRB 990308. 
The faintness of these hosts suggests that, regardless of 
the host luminosity, GRBs seem to be associated with star 
formation (SF). This correlation may allow GRBs to be 
used as a powerful tracer of star formation, provided the 
link between GRBs and SF is correct. It also implies that 
a significant amount of star formation is located in (op- 
tically) very faint galaxies. Their faintness may be due 
to dust extinction or a very steep faint-end slope of the 
galaxy luminosity function. The SF/GRB correlation may 
therefore be the basis of a unique and new way of finding 
star-forming galaxies at high redshifts (independent of the 
otherwise unavoidable surface-brightness bias). 

Bloom et al. (2002) present a comprehensive study of 
20 GRB hosts, including the three reported here, based 
on the same HST data and mainly Keck imaging. For 
GRB 980329 they found an OA position more centered 
on the main host complex compared to our localisation. 
For GRB 980519 their OA localisation error is similar 
to ours but they do not identify the northern object as 
the host center. The reason for this is not entirely clear, 
but from their thumbnail image (their Fig. 2) it seems 
the northern object is much less significant than in our 
image and likewise the northern most blob/knot, which 
is prac tically coincident w ith the OA location. For GRB 
990308 [Bloom et al. (2002] ) obtain an OA localisation with 
significantly larger errors than ours, and identify a large 
galaxy in the outskirts of their error circle (indicated in 
Fig. ||). Our host/OA localisation corresponds to the ob- 
ject^) seen left of the center of the error circle in their 
image. In summary, we identify different host candidates 
in two out of the three faint hosts investigated. These 
three hosts are also among the four most extreme outliers 
(the ones with the largest OA to host-center offset) in the 
Bloom et al. sample. Adopting our host identifications and 
OA localisations, the sample does not have obvious out- 
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liers. This is quantified by re-computing P(n c h anC c = 0), 
defined in Bloom ct al. (2002), representing the probabil- 
ity that none of the host identifications of the Bloom et 
al. sample are randomn galaxies (unrelated to the GRB). 
Bloom ct al. (2002[) found P(n chancc = 0) = 0.483, but 



by using our OA localisations and host identifications we 
find P(n c hancc = 0) = 0.78. This means that by using our 
results it is unlikely that any of the 20 hosts in the Bloom 
et al. sample are false identifications. 
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